Monte-Carlo method for multiple parametric integrals 
calculation and solving of linear integral Predholm equations 
of a second kind, with confidence regions in uniform norm. 

E.Ostrovsky a , L.Sirota b 

a Corresponding Author. Department of Mathematics and computer science, 
Bar-Ilan University, 84105, Ramat Gan, Israel. 

E - mail: eugostrovsky@list.ru 



b Department of Mathematics and computer science. Bar-Ilan University, 84105, 

Ramat Gan, Israel. 

E - mail: sirota@zahav.net.il 



Abstract. In this article we offer some modification of Monte-Carlo method for 
multiple parametric integral computation and solving of a linear integral Fredholm 
equation of a second kind (well posed problem). 

We prove that the rate of convergence of offered method is optimal under nat- 
ural conditions still in the uniform norm, and construct an asymptotical and non- 
asymptotical confidence region, again in the uniform norm. 

Key words and phrases: Kernel, Linear integral Fredholm equation of a sec- 
ond kind, Monte-Carlo method, random variables, natural distance, Central Limit 
Theorem in the space of continuous functions, metric entropy, ordinary and Grand 
Lebesgue spaces, uniform norm, spectral radius, multiplicative inequality, confidence 
region, Kroneker's degree of integral operator, variance, associate, dual and conjugate 
space, linear functional, Dependent Trial Method. 

Mathematics Subject Classification 2000. Primary 42Bxx, 4202, 68-01, 62-G05, 
90-B99, 68Q01, 68R01; Secondary 28A78, 42B08, 68Q15, 68W20. 



1 Introduction. Notations. Problem Statement. 
Assumptions. 

We consider a linear integral Fredholm's equation of a second kind 

V(t)=f(t)+ I K(t,s)y(s) fi(ds) = f(t)+S[y}(t). (1.1) 
Jt 

Here (T = {t}, fi), s G T be me measurable space with a probabilistic: fi(T) = 1 
non-trivial measure /x, S[y](t) is a linear integral operator (kernel operator) with 
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the Immeasurable kernel K(-, •) : 

S[y](t) d = f J T K(t,s)y(s)^ds). (1.2) 

For example, the set T may be subset of the whole Euclidean space R d with non-zero 
Lebesgue measure fj,(ds) = ds. 

The case when the domain T dependent on the variable t may be reduce in general 
case after some substitution to the case of equation (1.1); it may be implemented, 
e.g., for the Volterra equation [17]. 

The equations of a view (1.1) appears in many physical problems (transfer equa- 
tion, potential theory etc.), in the reliability theory (renewal equation), in the nu- 
merical analysis, for instance, for computation of eigenfunctions and eigenvalues for 
integral operators etc. 

We denote as ordinary for arbitrary measurable function g : T — >■ R 

\g\ P = 



\g(s)\ p »(ds) 



v e [i,oo], 



g e L(p) \g\ p < oo, 
where we will write in the case p = oo 

\g\oo = vraisup \g(t)\. 
teT 

We define analogously for the random variable £ 

|£| P =[E|£H 1/P ,p>l. 
Let us define the following important function, presumed to be finite \i — almost 
everywhere: 

R{s) = vraisup \K(t,s)\ (1.3) 
teT 

and we introduce also the so-called natural distance, more exactly, semi- distance, 
d = d(t, s) on the space T : 

uj. \ de f ■ \K(t,x) - K(s,x)\ . . 

d(t, s) = J vraisup 1 v n , 1.4 

x eT R{x) 

so that 



\K(t,x) -K(s,x)\ < R(x) d(t,s). (1.5) 

We assume that the metric space (T, d) is compact set and the measure /i is 
Borelian. We suppose also the function f(t) is d — continuous / : T — > R 1 . 

It follows from the inequality (1.5) that the function t — > K(t,x) is also d — 
continuous for /i — almost everywhere values x; x G T. 
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For the integral operator S[-] the m th , m > 2 power of S : S m [-] may be calculated 
as usually 

f II — ']_ Til, 

S m [g](t) = f m K(t, Sl ) J] K( 8j ,8 S+1 ) g(s m ) H^dsi). (1.6) 

Tm j=i i=i 

We will consider the source equation in the space C (T) = C (T, d) of all d — 
continuous numerical functions g : T — >■ i? with ordinary norm 

Ibll = sup|y(i)| = bloo- (1.7). 
teT 

Recall that the norm of linear integral operator S[-} (1.2) in this space may be 
calculated by the formula 

||5||^ sup ||%]||/||^|| = sup / \K(t,8)\n(d8). (1.8) 

A spectral radius r = r(L) of a bounded linear operator L : C{T) — > C(T) is 
defined by the expression 

r(L) = lmi m _ >00 II L m II 1/m = lim ||L m || 1/m . (1.9) 

m— s-oo 

More detail information about spectral radius of operator see in the classical books 
of N.Dunford, J.Schwartz [9], chapter VII, section 3; [10], chapter IX, section 1.8. 
For instance, it is known that r(L) < \ \L\\. 

We define for the kernel operator S[-] of a view (1.2) the so-called Kroneker's 
power U = as an kernel integral operator by the following way: 

U[g](t) = J T K 2 (t,s) g(s) ^(ds). (1.10) 
We assume (the essential condition!) that 

p d M r {U)=r(sW)<l. (1.11) 
Note that in the considered case (1.11) 

Pi = f r(S)< Vp<1- (1-12) 

Example 1. Let the set T be closed interval T = [0, 1] equipped with classical 
Lebesgue measure. Suppose the function (t, s) — > K(t, s) is continuous and denote 

7 = 1 (K) = sup \K(t,s)\. (1.13) 

t,sGT 

If 7(fT) < 1, then evidently p^S) < 7 < l,p(S) < 7 2 < 1. 

Example 2. Let us consider instead operator S the Volterra's operator of a 
view 

V[<7](f)= r^(t,s) g(s) ds, 
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where T = [0, 1] and we suppose that the function (t, s) — > K(t, s) is continuous and 
we will use again the notations (1.13). Then 

\\V m \\ < — . 
ml 

Therefore 

Pl (V)=p(V) = 0. 

The (linear) operators V with the property r(V) = are called quasinilpotent. 

Notice that the Monte-Carlo method with optimal rate of convergence for linear 
Volterra's equation 

y(t)=f(t)+ f K(t,s) y(s) ds 
Jo 

is described, with an applications in the reliability theory, in an article [17]. We 
represent in this article some generalization of results of [17]; but we note that 
in this article is considered also the case of discontinuous function /(•,•) and are 
described applications in the reliability theory (periodical checking, prophylaxis). 
Example 3. We denote for arbitrary linear operator L 

r m (L)^\\L m \\. (1.14) 

If 

3/5i e (0, 1), Ci, Ai = const < oo =>• r m (S) < Cim Al /3 1 m , m = 1, 2, . . . , (1.15) 
then obviously p\ < Analogously, if 

3/3 G (0, 1), C,A = const < oo => r m (U) < Cm A f3 m , m = 1, 2, . . . , (1.16) 
then p < j3. 

We intend to prove that under formulated conditions that there ex- 
ists a Monte-Carlo method for solving of equation (1.2) with the classi- 
cal speed of convergence l/^/n, where n denotes the common number of 
elapsed random variables. 

Moreover, at the same result is true when the convergence is under- 
studied in the uniform norm. 

The letter C, with or without subscript, denotes a finite positive non essential 
constants, not necessarily the same at each appearance. 

The papier is organized as follows. In the next section we describe the numerical 
Monte-Carlo method for solving of integral equation (1.1) and prove the optimality 
of it convergence in each fixed point to; t e T. 
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In the third section we recall for reader convenience some used facts about Grand 
Lebesgue Spaces of random variables and random processes and obtain some new 
used results. 

Fourth section is devoted to the Monte-Carlo computation of multiple parametric 
integrals. The fifth section contains the main result of offered article: confidence 
region for solution of Fredholm's integral equation in the uniform norm under the 
classical normalizing. In the next section we built the non-asymptotical confidence 
domain for multiple parametric integrals and for solution of Fredholm's integral 
equation in the uniform norm. In the 7 th section we consider some examples to 
show the convenience of using of offered algorithms. 

In the next section we offer the Monte-Carlo method for derivative computation 
for solution of Fredholm's integral equation again in the uniform norm with optimal 
rate of convergence. 

The 9 th section included some additional remarks. The last section contains 
some results about necessity of conditions of our theorems. 



2 Numerical method. Speed of convergence. 

A. Deterministic part of an error. 

Let e > be arbitrary "small" positive number. The solution of the equation 
(1.1) may be written by means of the Newman's series 

oo 

v(o = /(o + E<n/](*)- (2-i) 

m=l 

Let e G (0, 1/2) be a fix "small" number. We introduce as an approximation for 
the solution y(-) a finite sum 

y {N) (t) = £ s m [f](t), (2.2) 

m=l 

where the amount of summands N = N(e) may be determined from the condition 

oo 

E \\S m [f]\\<e. (2.3) 

m=N+l 

If for example the operator S[-] satisfies the condition (1.15), then the number 
N = N(e) may be find as the minimal integer solution of the inequality 



Ci\\f\\ £™ Al /3r < e (2-4) 

JV+1 



in the domain 



iV > argmax m Al /3™ 



log A 
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The asymptotical as e — > 0+ solution N ~ N = N (e) of the equation (2.4) has a 
view 



N = log(Ci|log/3i|/ei) + Ax [log(C 1 | \o g p 1 \/e 1 )/\ log AO , (2.5) 
ei = e/ (Ci 



B. Probabilistic part of an error. 

We offer for the calculation of the value y( N \t) at the fixed point t E T the 
Monte-Carlo method. Let us introduce a following notations. 

oo 

R a = R a {U) = ^2 k a rl^ 2 {U), a = const G (— oo, oo); 

k=l 
N 

R a {N, U) = J2 k a rl /2 (U), a = const e (-oo, oo). 
k=i 

Further, for the m— tuple x — (x±, X2, ■ ■ ■ , x m ), m = 2, 3, ... we define 

i?W [f](t, x) = K(t, x 1 )K(x 1 ,x 2 ) . . . K(x m . u x m )f{x m ). (2.6) 
We define in the case m = 1 

K^[f}(t,x) = K(t,x 1 )f(x 1 ). (2.7) 
Let ^ be a double sequence independent \x distributed random variables: 

P(g ) eA)=tL(A), 

where A is arbitrary Borelian subset of the space T. We introduce the following 
random vector for the integer values m — 1, , 2, . . . : 

pU) _ sti pi pi \ 

Sm LSI) S2) • • • i >>mi " 

Let us denote 

N 

h = (n(l), n{2), . . . , n(iV)), £(ra) d = ^ m • n{m), 

m=l 

where N = N(e), n be any N— tuple of positive integer numbers: n(j) = 1,2, ... . 
We consider the following Monte-Carlo approximation for the multiple integral 
S m [f] : 

1 n ^ -» 
Sn U) [f] = S™ {j) [f](t) = £ K^[f](t,^) (2.8) 
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and we offer correspondingly the following approximation y^\t) for the variable 

y {N \t): 

N 

y { "\t) :=/(*) + £ ^)[/K0- 

i=i 



(2.9) 



Note that the approximation (2.8) for the (multiple) parametric integrals was intro- 
duced by Frolov A.S.and Tchentzov N.N., see [14], and was named "Dependent Trial 
Method". It was proved under some hard conditions that the rate of convergence 
of this approximation in the space of continuous functions is optimal, i.e. coincides 
with the expression 1 / y/n, where n denotes the amount of all used random variables. 

Theorem 1.1 We assert under formulated conditions that for arbitrary e G (0, 1) 
and for all sufficiently great values n; n — > oo there exists the tuple h = n(n) for 
which 



and 



sup min Varjyf (e) \t) -y {N \t)) < 

t£T n:B(h)<n ^ ' 



2 R 1/2 (U) R- 1/2 (U) 



1 | CiiP) 



n 



(2.10) 



sup min Var (y^^\t) — 

t£T n:B(h)<n ^ 



> 



R 1/2 (U) R- 1/2 (U) 



1 C 2 (U) 
n n 2 



where the positive finite constants C\ = C\(U), C 2 = C 2 (U) does not depend on the 
n, e and /. 

Remark 2.1. Note that under condition p(U) the values R 1 / 2 (U) and R_i/ 2 (U) 
are finite. 

Proof of the theorem 1. Without loss of generality we can suppose ||/|| = 1. 
We conclude that the variance of each summand in (2.9), i.e. the expression 

v(m,n( 3 m) :=Var (%[/](*)) 
may be estimated as follows: 



Therefore, 



supf (m, n(j))(t) < r m /n(m), m — 1, 2, . . . 



N r- 



t£T 



(2.11) 



On the other hand, the common amount of used random variables in the 
formula (2.9) is equal to the expression 
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N 



B{n) = m ' n(m). 



(2.12) 



m=l 



Let us consider the following constrained extremal problem: 

$(n) — > min / 1 B(n) = n. 

We obtain using the Lagrange's factors method the optimal value no = 
{n (l), no(2), . . . , } of the tuple h for the problem (2.13) has a view: 



(2.13) 



n (m) 



nrU 2 (U) def 



Ru 2 (N,U) ^ 



6{m) n, 



up to rounding to the integer number, for example, 

nrU 2 (U) 



n (m) = 1 + Ent 



R 1/2 (N,U) 
1 + Ent [6{m) n] , 
where Ent (a) denotes the integer part of the positive number a and 



(2.14) 



6{m) = 9(m,N,U) 



R 1/2 (N,U) ^ 



The minimal value of the functional <£>(n) under condition B(n) = n may be 
estimated as 



min$(n)/[B(ra) = n] < R 1/2 (N, U) R- 1/2 (N, U) 



n 



rr 



< 



Ri/2(U) R- 1/2 (U) 



n 



The lower bound provided analogously. 
This completes the proof of theorem 1.1. 

Example 2.1. Assume that the sequence rk(S) satisfies the condition (1.15). 
As long as as x — > 1 — and (5 = const > —1 



k=i 



2^ +1 r(/3 + 1) 

llogrrl^ 1 ' 
we have as a = const > — 1/2, Pi — >■ 1 — 

supVar( y f^W- y W(t)), 



Ci 



t<=T 



2 2 2Q+2 r(a + 3/2)r(a + l/2) 



log A 



\2a+2 



1 , Ci(U) 



n 
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Remark 2.1. Note that in the case T C R d for each random number £ gen- 
eration are used in general case d uniform distributed in the interval [0, 1] random 
variables. See, e.g. [7], [17]. 

3 Banach spaces of random variables 

Pilcrow A. Banach spaces of random variables with exponentially decreasing tails 
of distributions. ("Exponential" level). 

In order to formulate our results, we need to introduce some addition notations 
and conditions. Let = 0(A), A G (— Ao,Ao), Ao = const e (0, oo] be some even 
strong convex which takes positive values for positive arguments twice continuous 
differentiable function, such that 

0(0) = 0, // (O) > 0, lim 0(A)/A = oo. 

A— >Xq 

We denote the set of all these function as $; $ = {0(-)}- 

We say that the centered random variable (r.v) f = belongs to the space 
5(0), if there exists some non-negative constant r > such that 

VAe(-A ,A ) => Eexp(A0 < exp[0(A r)]. (3.1). 
The minimal value r satisfying (4) is called a -8(0) norm of the variable £, write 

M\B{4>) = inf{r, r > : VA Eexp(A0 < exp(0(A r))}. (3.2) 

For instance, if 0(A) d = 02(A) = 0.5A 2 , A e R the space B{4> 2 ) is called subgaussian 
space and is denoting ordinary -6(02) = sub = sub(f2) in accordance to Kahane [21]; 
the (centered) random variables from this space are called subgaussian. 

The norm in subgaussian space of a random variable £ will denoted ||£|| sub. 

The important example of subgaussian random variables (r.v.) are centered 
Gaussian (normal) variables; indeed, if r.v. Law(£) = iV(0,<7 2 ),o- > 0, then 

||£||sub = (7. 

If a centered r.v. £ is bounded, then it is also subgaussian and 

Hfllsub < |f |oo := vraisup |f|. 
For instance, the Rademacher's r.v. £ : 

p(e = i) = p(e = -i) = 1/2 

is also subgaussian and ||f || sub = 1. 

It is proved in the article [4] that the space sub(fi) is Banach space. The centered 
random variable £ belongs to the space sub(f2) and has a norm r = ||f || sub, r > 
if and only if 

VA e R Eexp(Af) < exp (o.5AV) . 
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More details about the space sub(fi) see in the article [4]. 

The spaces B(<f)) are rearrangement invariant in the terminology of a book [1], 
chapter 2 and 3; [25], chapter 3; are very convenient for the investigation of the r.v. 
having a exponential decreasing tail of distribution, for instance, for investigation of 
the limit theorem, the exponential bounds of distribution for sums of random vari- 
ables, theory of martingales, non-parametrical statistics, non-asymptotical proper- 
ties, problem of continuous of random fields, study of Central Limit Theorem in the 
Banach space etc., see [3], [23], [55], [56], [57], [58], [30], [31], [33], [34], [35], [39], [40], 
[41],[42],[43], [55]. 

The generalization of this spaces on the case = oo is considered, e.g., in 

[6], [12], [13], [18], [19], [26], [44], [45] etc. 

The space B((f)) with respect to the norm || • and ordinary operations is 

a Banach space which is isomorphic to the subspace consisted on all the centered 
variables of Orlichs space F, P), N(-) with N — function 

N(u) = exp(0*(w)) - 1, 

where 

4>*{u) = / sup(Aw-0(A)). 

A 

The transform — > 0* is called Young-Fenchel transform. The proof of consid- 
ered assertion used the properties of saddle-point method and theorem of Fenchel- 
Moraux: 

see [24], chapter 1. 

Let £ be centered random variable such that its moment generating function 

A Eexp(A£) 

is finite in some neighborhood of origin: |A| < A , A = const, < A < oo. 

Here A may be complex; in this case the moment generating function is analytical 
inside the circle |A| < A . 

The finiteness of moment generating function is equivalent the following moment 
inequality: 

\Z\ P <Cp, P >1. 

The natural function = 0^(A) for the variable £ may be introduced by a formula 

€ (A) = logEexp(A0. 

It is obvious that 0^(-) G 

Analogously is defined a so-called co-transform v — > v* : 

v*(x) = inf (xy + v(y)). 
i/e(o,i) 
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The next facts about the -8(0) spaces are proved, for instance, [23], [30], p. 19-40: 
1. f G B{(j>) <=> E£ = 0, and 3C = const > 0, 



U(£,x) < exp(-(f)*(Cx)),x > 0, 
where U (£, x) denotes in this article the tail of distribution of the r.v. £ : 

C/(£,x)=max(P(£>x), P(£<-a;)), x > 0, 

and this estimation is in general case asymptotically exact. 

Here and further C,Cj,C(i) will denote the non-essentially positive finite "con- 
structive" constants. 

More exactly, if A = oo, then the following implication holds: 

lim 0- 1 (logEexp(AO)/A = K G (0, oo) 

A— s-oo 

if and only if 

lim { ( j ) *)- l (\\ogU^,x)\)/x=l/K. 



Here and further denotes the inverse function to the function / on the left-side 

half-line (C, oo). 

Let £ = = €(t,u) be some centered random field, t G T. The function 
0(A) = 0£ (A) may be constructive introduced by the formula 

0(A) = O (A) = f logsu P Eexp(A£(t)), (3.3) 

ter 

if obviously the family of the centered r.v. {£(£), t G T} satisfies the uniform 
Kramers condition: 

G (0, oo), sup U (£(t), x) < exp(— /i x), x > 0. 

In this case, i.e. in the case the choice the function 0(-) by the formula (3.3), we 
will call the function 0(A) = 0o(A) a natural function for the random process £(£). 

Pilcrow B. Grand Lebesgue Spaces of random Variables. ("Power" level.) 

Let ip = i/j(p) be function defined on some semi- open interval of a view 

1 < p < 6, 

where obviously b = const, 1 < b < oo, is continuous and is bounded from below: 
inf pG ( l b ) ip(p) > function such that the function 

w(p) = w^(p) = plogip(p),p G (1,6) 

is downward convex. We will denote the set of all such a functions as \I/ = \I/(1, 6) : 

if) G 
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We define as ordinary 

supp ip = [1, 6) 

and put for the values p > b — > if)(p) — +00 in the case b < 00. 

It may be considered analogously the case of closed interval p G [1,6], 6 G (1, 00]; 
then supp-0 = [1, b] and for all the values p > b — > ip{p) = +00. 

We introduce a new norm (the so-called moment norm) on the set of r.v. defined 
in our probability space by the following way: the space G(ip) consist, by definition, 
on all the centered r.v. with finite norm 

m\G^) d ^ f sup la/iM, ia = El/ ^r- (3.4) 

Remark 3.1. Note that in the case supp-?/' = [1,6], 6 = const G (l,oo) the norm 
1 1£| | G(ip) coincides with the ordinary Lebesgue norm |£| 6 : 

Indeed, the inequality < ||£||Gr(^) is evident; the inverse inequality follows from 
the Lyapunov's inequality. 

Let us introduce the function 

x(p) = XAP) = ^T^-y P>1- 

It is proved, e.g. in [23], [30], chapter 1, section (1.8) that the spaces -8(0) 
and G(x) coincides: 5(0) = G(x) (set equality) and both the norm || • \ \B((f)) and 
|| ■ \\G(x) are equivalent: 3d = Ci(0),C 2 = C 2 (0) = const G (0,oo), V£ G B(<j>) 

m\G(x)<c 1 m\B^)<c 2 m\G( X ). 

Conversely, for arbitrary function ifj(-) G \I/(1, oo) for which 

_ log^(p) 

hmp^oo— < 1 

logp 

may be defined the correspondent function = 0(A) as follows: 

, A > A = const > 0. 

Recall that at |A| < A 0^(A) x C\ 2 . 

The definition (3.4) is correct still for the non-centered random variables £. If for 
some non-zero r.v. £ we have ||£||G(-0) < oo, then for all positive values u 

P(|£| > u) < 2 exp (- W ;(«/(C73 ||e||GW))) • (3-5) 
and conversely if a r.v. £ satisfies (3.5), then ||£||G(-0) < oo. 



A 
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The definition (3.4) is more general as (3.1). Indeed, if a r.v. £ belong to some 
space B(4>), G then Vp G (l,oo) |£| p < oo. The inverse inclusion is not true, 
e.g., for the symmetrical distributed random variable £ for which 

P(|C| > u) = exp (-m a ) , u > 0, (3.6) 

where A = const G (0, 1). 

Further, let £ be any r.v. such that for some b = const > < oo. The natural 
choice of the function ^(p) f° r the r.v.£ may be defined by the formula 

Mp) = \Z\p, P--\Z\ p <oo. (3.7) 

Remark 3.2. Note that: 

A. The r.v. £ is bounded if and only if 

linip^oo^p) < oo. 

B. The r.v. £ satisfies the Kramer's condition if and only if 

lim^oolog^(p)/logp < 1. 

C. The r.v. £ obeys all the exponential moments, i.e. 

VA G R Eexp(A£) < oo 

if and only if 

linip^log^pVlogp < 1. 
For instance, for the r.v. ( in (3.6) the natural function ip^(p) has a view 

Mp) = ICIp = 2 1/p r 1 ^( p /A + 1), i < P < oo. 



Note that as p — > oo 



^c(p) ~ (p/(Ae) 1/A - 



Pilcrow C. Non- asymptotical bounds of distributions in the classical CLT. 
Let us define for all the functions G $ 

0(A) = sup n(f)(X/y/n). 

n=l,2,- 

For example, let 

0(A) = r (A), |A| > 1 r (A) = Ci|A| r , r = const > 1; 

then 

0r(A) X max (r,2)(A), |A| > I. 
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We denote also for all the functions ip £ ^(2, b) 

^(p) = Cq 1 -p i/j(p)/ log p. 

The provenance, calculation and exact value of the constant C ~ 1.77638... is 
described in [45]. 

Obviously, if b < oo, then ip{p) x ip{p). It is not true in the case when b = oo. 
The probabilistic sense of introduced functions is following. Let £ G B(<p) and 
let rj G Er^ = 0. Let also = 1,2,... be independent copies of £ and let 

be independent copies of rj. Then 



Bup||n- 1/2 E^)IW)<||eiW), (3.8a), 



i=i 



(exponential level), 



supWn-^^vUm^) < \\v\\G(i>). (3.86) 

3=1 



(power level). 

Remark 3.3. It is important to notice that the exponential bounds for the tail 
behavior for the sums of independent random variables (3.8a) does not be obtained 
from the power estimations (3.8b) and conversely the power estimations (3.8b) does 
not be obtained from exponential estimations (3.8a). Let us consider the following 
examples, see [30], p. 55-57. 

Let i — 1,2, ... be a sequence identical distributed centered r.v. with the 

following tail function: 

U(£i,u) = exp (— u r ) , r = const > 0, u > 2. 

Denote 



P r {u) = sup P n 



-1/2 



i=i 



> u , u > 2. 



From the relation (3.8a) follows the estimation 



P r {u) < exp (-C(r) M min(r ' 2) ) , u > 2, (3.8c) 

only under the condition r > 1. 

But from the relation (3.8b) follows the inequality 

P r (u) < exp (-Ci(rK /(r+1) [log M ] C2(r) ) , u > 2. (3.8d) 

Note that the inequality (3.8c) is more exact that (3.8d), but only in the case r > 1. 

In the case r G (0, 1). the r.v. does not belong to the any B(<p) space, G $; 
therefore the relation (3.8d) has advantage in the considered variant. 
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Pilcrow D. Continuity and compactness of random fields. 

Let £(t) = £(t,u), t G T be a separable centered random field. The "con- 
structive" introduction of the ip = ip^(p) function for the random field £(•) may be 
follows: 

^(p) = sup |£(f) | p , (3.9) 

if it is finite for some p > 1. 

The natural function ^(A) for the £(•) is defined as follows: 

0c(A) = maxsuplogEexp(A fi £(£)), (3.10) 
if reasonably the last function is finite for some non-trivial interval 

A G (— A , A ), A = const > 0. 

Evidently, ip^(-) G ? (-) G $. 

Let us denote for arbitrary function tp G ^ 

tie/ 

u*(a:) = u*,^(a;) = inf (xy + log^(l/y)). (3.11) 
ye(o,i) 

M.Ledoux and M.Talagrand [55], chapter 2 write instead our function 
exp (— <f>*(x)) some Youngs function ty(x) and used as a rule a function *ff(x) = 
exp(— x 2 /2) (the so-called subgaussian case). 

Without loss of generality we can and will suppose 

sup[iie(*) iw)i = i, 

ter 

(this condition is satisfied automatically in the case of natural choosing of the func- 
tion 4> : 0(A) = 0o(A) ) and that the metric space (T,d) relatively the so-called 
natural distance (more exactly, semi-distance) 

= (3-12) 

is complete. 

Recall that the semi-distance d = d(t, s), s,t G T is, by definition, non-negative 
symmetrical numerical function, d(t,t) = 0, t G T, satisfying the triangle inequality, 
but the equality d(t, s) = does not means (in general case) that s — t. 

For example, if £(£) is a centered Gaussian field with covariation function 
D(t, s) = Ef (t) f (s), then O (A) = 0.5 A 2 , \ E R, and rf(t, s) = 

Mi) -£(s)||5(0o) = VVarf^t) - £(s)] = ^D(t, t) - 2D(t, a) + s). 

Let us introduce for any subset V, V C T the so-called entropy H(V, d, e) = 
e) as a logarithm of a minimal quantity N(V, d, e) = N(V, e) = iV of a balls 
5(\/,t,e), teV : 

S(V, t, e) d = {s, seV, d(s, t) < e}, 
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which cover the set V : 

TV = min{M : 3{U}, i = 1, 2, , M, U G V, V C U^V, t i? e)}, 
and we denote also 

B(ti,e) = logiV; S(t ,e) = f S(T,t ,e), H(d,e) = f H(T,d,e). 
A capacity M = M(V, d, e) of the set V is the maximal number of disjoint balls 
B(ti, e) = {s,se V, d(U, s) < e}, ti, s G V, i = 1, 2, . . . , M, 

B(t h e) n B(t,, e) = 0, i± j, i, j = 1, 2, . . . , M, 
subsets of the set V : 

uf =1 B(U,e)cV. 

Denote 

H(e) = H(T,d,e), L(e) = logM(T,d,e), 

M e ) = j n | ^( B (t, e )), M e ) = sup fj,(B(t,e)). 
It is known [30], p. 95-97; [60], p. 5-9 that 

H(2e) < L{e) < H{e). 
The next fact allows to estimate the values H(e),L(e) : 

"< T > <M(T, C )<^ 



h + (2e) ~ v ' ' ~ h-(e)' 

If for instance the set T is bounded open subset of the space R d equipped with a 
distance d(t, s) for which 

d\t - s\ a < d(t, s) < C 2 \t - s\ a , a = const G (0, 1], 



\t - s\ 



\ J2( f k ~ S k) 2 , 
\ k=l 



then 

N(T,d,e)-C,(T,a) e~ d ' a , e G (0, 1). 

It follows from Hausdorf 's theorem that Ve > =>- H(V, d, e) < oo iff the metric 
space (V, d) is precompact set, i.e. is the bounded set with compact closure. 

Let £(£), £ G T be a separable numerical random field such that for some function 
ip G \& sup tgT ||£(£)||Cr(^) = 1. We introduce so-called natural distance 
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d4t,s) = \m-m\\GW. (3.13) 

Denote 

v(y) = log^(l/y), v*(x) = inf + 

y:ip(l/y)<oo 

Lemma 3.1. a If the following integral converges: 

K(a<b) ■= I x<t>( H (T,dd„x)) dx < oo, 
Jo 

where 

°<f> = sup | \£(t)\\B(<f)) < oo, 
teT 

then (see [23], [30], chapter 3, section 3.4) 

P(t(-)eC(T,d*)) = l 

and for the values 

u>2K ((T ) 

we aave 

F sup o$)u ) ) . (3.14a) 

Lemma 3.1.b If the following integral converges: 

[ exp (w*(log(2if(T, gL,x)))) dx < oo, 
Jo 

then (see [30], chapter 4, section 4-1) 

P(£(-)eC(T,aV)) = l. 



Moreover, 



sup \m)<z(8), 

t,s:d^(t,s)<S 



where 

Z(5) = Z{i), 5) = f 9 / exp (u„(log(2if (T, cL, a;)))) dx. (3.146) 

JO 

As a slight consequence: let us denote 

<ty = sup||£(t)||G(^) < oo. 

ter 

We have for arbitrary fixed non-random value to & T using triangle inequality and 
taking in (3.14b) 5 = : 
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sup|£(t)| < |£(*o)|+sup|£(f)-£(* )|; 
ter ter 

sup 1 \\G(iP) < \\ttto)\\G(ij) + Z M < 
teT 

<j^ + Z (cry,) d = Z = Z(ip) = Z(ip, e); 



P [sup \£(t)\ > u J < exp [-w], (u/Z{i)))] , u > 2Z(iP). (3.14c) 



teT 



For instance, let ip(p) — pP , f3 — const > 0. The condition (3.14b) may be written 



as 



/ H 13 (T, djf,, x)dx < oo. 

J 

Moreover, if 

ij){p) = p^ 1 log^ 2 (l +p), p> 1, ft = const > 0, 
then the condition (3.14) has a view 



f 1 H^(T,d^,x) log? 2 (1 + H(T,d^,x)) dx 

J 



< oo. 



Let now supp?/> = [1,6], 6 = const > 1. The condition (3.14) is equivalent to the 
famous Pizier's condition [37] 



f 1 N^ b (T,cL,x) dx < oo. 
Jo 



Further, let £ a (t), t e T, a e Abe arbitrary family of random fields, where A 
is arbitrary set, such that such that for some function ip e ^ 

supsup||a(t)||G'(^) = l. 
We introduce again the so-called natural distance induced by the family {£<*(•)} 

d4t,s)= f S up\\Ut)-Us)\\GW. (3.15) 

Denote 

v(y) = log^(l/y), u*(a;) = inf (xy + 

y:il>(l/y)<oo 

Lemma 3.2. Assume that for some t e T £/ie one- dimensional family of random 
variables {£ a (^o)} *s stochastically bounded: 

lim supP(|^(t )| > u) = 0. (3.16) 
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If the following integral converges: 



f 

Jo 



exp (u,(log(2if (T, cL, x)))) dx < oo, (3.17) 

10 

then the family of distributions /i a (-) generated by the random fields in the 

space C (T, d) : 

H*(B) = P (&(•) G 5) 

zs weakly compact. 

The detail explanation of the theory of weak convergence for probabilistic mea- 
sures in the metric spaces see in the monographs [32], [2], [49]. A main conclu- 
sion of this theory for the continuous random processes may be formulated as fol- 
lows. If the finite-dimensional distributions of the sequence of a random processes 
7] n (t), n = 1,2, .. . converge as n — > oo to the finite-dimensional distributions of a 
random process r)(t) and the distributions of the random processes 

fi n (B) = P(77„(-) G B) 

are weakly compact, then the distributions of arbitrary continuous functional F : 
C(T) — > R at the points rj n (-) converge to the distribution F(rj) : VA G R 

\h^~E exp(i\F(r) n )) = Eexp(i\F(i])). 

In particular: 

lim P(sup |77„(t)| >u) = P(sup \r](t)\ > u) (3.18) 
for all the points u, u G R in which the function 

u — >■ P(sup \r)(t)\ > u) 

is continuous. 

Pilcrow E. Central Limit Theorem ( CLT) in Banach space. 

Let (Q, A, P) be a probabilistic space and let B with a norm | ] | • 1 1 1 be a Banach 
space. We will say as ordinary that the (centered) random variable (, ( : Q — > B 
satisfies the CLT in this space, if the sequence 

n 

C„ = n" 1 / 2 ECi (3-19) 
i=i 

converges weakly in distribution as n — > oo to non-trivial Gaussian centered random 
variable Coo : 

hm Law(C n ) = Law(Cj. (3.20) 
Here the variables {Q} are independent copies (. 
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It is evident that the variable Coo has at the same covariation operator as the 
variable (. If (3.20) there holds, then for all positive values u 

JmP(|||C B |||>«) = P(|||Coolll>«)- (3-21) 

As well as in the one-dimensional case, the CLT in the Banach space will be used 
in the Monte-Carlo method for building of confidence region for estimated function 
in the Banach space norm. 

Pilcrow F. Multiplicative inequalities in Grand Lebesgue spaces. 

Let £ G G(i^i),r] G G(ip2), t — £ • T). We will study in this subsection the 
inequalities of a view 

\\r\\GM = \m\GM < C\\£\\GM \\v\\GM. (3.22) 
("power" level) or analogously 

IM|B(&) = 11^115(03) < CWtWBfa) \\r,\\G(<h). (3.23) 

("exponential" level). 

It is convenient to continue arbitrary function ip — i]){p) as follows: 

Vp ^ supp(^) ip(p) = +oo. (3.24) 

A. Independent case, power level. 

Note first of all that if the r.v. £,77 are independent, then 

\Cv\ P = \C\ P Hp- 

Therefore, if £ G G(V'i), V ^ G(ip 2 ), then 
Proposition 3. A. 

imm^.^^miGM-MiGM. (3.25) 

B. Dependent case, power level. 

We do not suppose in this pilcrow the r.v. £,77 to be independent. Let again 
£ G G(ipi), r] G G{jp2)- We introduce the following operation for two functions 
V'i(')) ^2(0 from the set \I/ : 

A o Mr) = infill ipr) ■ Mrp/(p ~ I))}- (3.26) 

p>i 

If supp^i = [1, 61), supp-^2 = [1,&2), then 

supp 0^2(0) = t 1 '^), 

where 

63 = 6163/(61 + 62). (3.27). 
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Note that if (61 - l)(b 2 - 1) > 1, then 63 > 2 and if (61 -2)(6 2 — 2) > 4, then 63 > 2. 
It is known [26] that in the general, i.e. dependent case 
Proposition 3.B. 

||^l|G(^o^)<||e||G(^i).|H|G(^). (3.28) 
C. Independent case, exponential level. 

Let £ G B(4>i) and r\ G B((f) 2 ) be independent r.v. and 01,2(0 be two functions 
from the set $ such that 

^Wt^Wj < 2 (3 . 29) 
logp 

We define a new operation (commutative and associative) 3 (p) = 0i 02 (p) for 
the functions 0i(O and 2 (O as follows: 

3 (p) = [0i~» 2 '(P)/^ 1 , P > 1- (3.30) 

Recall that at |A| < 1 =4> 3 (A) = C • A 2 , where the constant C must be choose 
such that the function A — > 03(A) is continuous. 
Note that by virtue of independence 

it l / ^2 P P ^ n P 

01 ip) 02 \P) 03 (P) 

therefore under condition (3.29) 
Proposition 3.C. 

||£77l|5(V3)<C 4 ||£||S(0i)-|H|S(02). (3.31) 

As a consequence: if the r.v. £1, £2, • • • , £fc are mutually independent and £j G B((f>j), 
then 



116 6 ...6||5(0 (fe) )<Qll6l|5(0i) ||6|| J B(02)...||6|| J B(0 fc ), (3.32) 

where 

{fc) = (((0i0 0i)0 2 )...0 fe ). (3.33) 
If for example the r.v. £, 77 are subgaussian and independent: 

1 1£| I sub = o"i, I |r/| I sub = a 2 , &i < 00, 

then 

|£| P < C 4 <7i ^/p, \rj\ p < C 4 a 2 y/p, 
and the r.v. r = £ • rj satisfies the following moment condition: 
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\r\ p < C 5 Gi a 2 p, p E [1, oo), 
or equally the r.v. r satisfies the Kramer's condition. 
D. Dependent case, exponential level. 

Let again £ e £>(0i) and 77 e -8(02) be arbitrary r.v. and 0i,2( - ) be two 
functions from the set $. We define a new operation (commutative and associative) 
4>a{p) — 0i <S> 02 (p) for the functions 0i(-) and 2 (-) as follows. Denote 

^■(p) = .(-f), r » i = 

0} >{p) 

*P 4 (p) = [fa fa](p), 

then we define 



i(p) = [0i®0 2 ](p) 



P 



^ 4 (p) 



(-1) 

(3.34) 



if obviously for some 

6 > 1 =>. (f) 4 (b) < 00. (3.35) 

Note that for the CLT in the Banach space we need to assume b > 2. 

The function 4 (p) = [0i <8>02](p) has the following sense. If the condition (3.35) 
is satisfied, then 

Proposition 3.D. 

||ei7lW4)<C7||e||S(^).|H|S(^). (3.36) 

As a consequence: if the r.v. £1,^2, • • • ,6 are mutually independent and 6 G 
B((j)j), then 



6 ...6l|5(0( fe) )<C 4 1l6l|5(0i) ||6I|5(0 2 )... 116115(0/0, (3.37) 

where 

0(fe) = (((0i ® 0i) <8> 2 ) • • • ® 0fc). (3.38) 

Let us return for instance to the subgaussian case. If the r.v. 6 V are subgaussian 
(and arbitrary dependent): 

||£|| sub = a 1 , \\t]\ \ sub = a 2 , (?i < 00, 
then the centered r.v. r = £ ■ 77 — E(£r?) satisfies the following moment condition: 

\r\ p < Cq cti (T 2 p, p e [1, 00), C 6 > C 5 , 
or again the r.v. r satisfies the Kramer's condition. 
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Let us consider more essential but more exotic example. Let {£} be symmetri- 
cally distributed r.v. with the following tail behavior: 

P(f > x) < exp (-Ci e x ) , x > 1, (3.39) 

and put 

k k 

v =ffc = n^ _E n^' 
i=i i=i 

where & are independent copies of £. It follows from proposition 3.D after some 
calculations 

U(v k , x) < exp (-C 2 (Jfe, d) e xl/k ) , x > 1, (3.40) 

In order to illustrate the inequality (3.40), let us consider the following example. 
Let rji — £, i — 1,2, ... , ( > 0, and 

P(C > x) = exp (- e x ) , x > 1, 
and let v = nf=i Vi = \C\ k ', then 

P(i/ > x) = P (|C| > x 1/fc ) = exp (-e^) , x > 1. 

Remark 3.D. Let us consider the -B (0 e ,«) , K > 1 space as a -8(0) space with 
the correspondent function 

e , K (A)x|A|.log K (2 + |A|),|A|>l, (3.41) 

or equally 

log[0:»] x (CWX 1 '") ,x>l, 
Let also £ = t G T be separable random field such that 

SU P ||e(t)|| J B(0 e , K ) = l. 

Introduce the following (finite) distance 

p e , K (t, s) = \m ~ m\\B (<f>e, K ) ■ (3.42) 
If the following integral converges: 

Ie, K = f log 2K (1 + H(T, p e , K , x)) dx < oo, (3.43) 

J 

then 

P(£(-)eC(T,p e , K )) = l (3.44) 
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and 

P(sup \£(t)\ >u)< exp (-C 3 exp \C 4 (I e , K )x 1/K ] ) , x > 1. (3.45) 
ter v L J/ 

Pilcrow G. Z>ua/ spaces for Grand Lebesgue spaces. 

We describe briefly in this pilcrow the dual (conjugate) spaces to the Grand 
Lebesgue Spaces. 

We do not need to suppose (only in this Pilcrow G!) the finiteness condition 
fj.(T) = 1; it is sufficient to entrust on the measure /i instead the condition of sigma- 
finiteness; and suppose also the triplet (T, is resonant in the terminology of 
the classical book [1]. This imply by definition that either the measure \i is diffuse: 

VAgS, < fi(A) < oo => 3B c A, fi(B) = fi(A)/2 

or the measure \i is purely discrete and each atoms have at the same (positive) 
weight. 

Further, let (a, b), 1 < a < b < oo be a maximal open subset of the support of 
some function ijj(-) G We will consider here only non-trivial case when 

max(^(a + 0),^(6-0)) = oo. (3.45) 

Note that the associate space for GLS spaces are describes in [48]; a less general 
case see in [12], [13]. 

Namely, let us introduce the space DG(vp) = DGtp(a, b) consisting on all the 
measurable functions {g}, g : T — >■ R with finite norm 

\\g\\DG(il)) = inf inf ^(p(k))\gk\p(k)/(p(k)-i)}, (3.46) 

where interior inf is calculated over all the sequences {gj}, finite or not, of a mea- 
surable functions such that 

3 

and exterior inf is calculated over all the sequences {p(k)}, belonging to the open 
interval (a, b). 

An action of a linear continuous functional l g , g G DG(ip) on the arbitrary 
function / G G(tp) may be described as ordinary by the formula 

l g (f) = J f(x) g(x) n{dx) 

with = II^IIDG^). 

Recall that the set of all such a functionals {l g } equipped with the norm \ \l g \ \ = 
\\g\\DG(ip) is said to be associate space to the space G(vp) and is denoted as usually 
G'(i:) - [G(i:)}>. 

Define the space G°(ip) as a (closed) subspace of a space G(ip) consisting on all 
the functions {g} from the space G(ip) satisfies the condition 
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lim -^r = 0, (3.47) 
and introduce correspondent quotient space 

GoW) = Gty)/G°ty). (3.48) 

It is known [48] that the space DG(tp) is associate to the space G(ip) and is dual to 
the space G°(^). 

It is easy to verify analogously to the case of Orlicz space (see [52], p. 119-121, 
[53], [28], [27] that each function / G G{jp) ma y be uniquely represented as a sum 

f = r + fo, (3.49) 

(direct sum), where f° G G°(ip), f G G (ip). 

More exactly, the function f a represented the class of equivalence under relation 

Therefore, arbitrary continuous linear functional L = L gv {f) on the space G(tp) 
may be uniquely represented as follows: 

L g , v (f) = J T f°(x) 9(x) V(dx) + f T f (x) u(dx), (3.50) 

where g G DG(ip), v G ba(T, E, //); 6a(T, E, //) denotes the set of all /imte additive 
set function with finite total variation: 

\u\ = \v\(T) = sup[i/(A) - i/(T \ A)] < co. (3.51) 

Age 

Note that the (generalized) measure v is singular relative the source measure (/,; 
therefore 

\\L g>v \\G*(i/>) = \\g\\DGty) + \v\. (3.52) 



4 Monte-Carlo method for the parametric inte- 
grals calculation 

We consider in this section the problem of Monte-Carlo approximation and con- 
struction of a confidence region in the uniform norm for the parametric integral of 
a view 

/(*)= / g(t,x)v(dx). (4.1) 

Here (X, E, v) is also a probabilistic space with normed: v(X) = 1 non-trivial 
measure v. 
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A so-called "Depending Trial Method" estimation for the integral (4.1) was in- 
troduced by Frolov A.S.and Tchentzov N.N., see [14]: 

n 

I n (t)=n- 1 J29(t,rii), (4-2) 
i=i 

where {r/i} is the sequence of v distributed: 

P( Vi e A) = u(A) 

independent random variables. 

We intend in this section to improve the result of the article [14] and its conse- 
quence, see [30], chapter 5, section 5.11. 

The modern methods of (pseudo)random variable generations are described, in 
particular, in [16]; see also [7]. 

We assume MteT g(t, •) e L 2 (X, v) : 

g 2 (t, x) i>(dx) < oo; 

then the integral I(t) there exists for all the values t;t e T and we can for any fixed 
point t eT use for an error evaluating the classical Central Limit Theorem: 

lim P(v^|J n (f ) - /(to) I < u) = $(u/a ) - $(-u/a ), (4.3) 
where as usually 

$(«) = (2tt)- 1/2 [ U exp(-z 2 /2) dz 

and 

a\t) = V a r(f(t,r ] ))= f g 2 (t,x)v(dx)-I 2 (t), 

o"o = cr(to). Let us consider now the problem of building confidence region for I(t) 
in the uniform norm, i.e. we investigate the probability 

P (n) (u) = f P(sup Vn\I n (t) - I(t)\ >u), u = const > 0. (4.4) 

On the other words, we use for construction of confidence region in the uniform norm 
in the parametrical case the Central Limit Theorem (CLT) in the space of continuous 
functions C(T) alike in the classical case of ordinary Monte-Carlo method is used 
customary CLT. 

Let us introduce as in the first section the following function, presumed to be 
finite v — almost everywhere: 

Q(x) — vraisup \g(t, x)\ (4.5) 
ter 

and we introduce also the so-called a new natural distance, more exactly, semi- 
distance, j3 — j3(t, s) on the space T : 
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at* \ def . \g(t,x) — g(s,x)\ ,, 
(5{t, s) = J vraisup n , \ , (4-6) 



so that 



\g(t,x)-g(s,x)\<Q(x) 0(t,s). (4.7) 
Let us consider the centered random process 

g°(t,x)=g(t,x)-I(t); Eg°(t,r))=0, 
and define the ip— functions as follows: 

HP) = MP) = sup \g%v)\ P = supE 1 /P|^°(i,77)| 1 ' = 

teT teT 



I \g(t,x)-I(t)\ p u(dx) 



1/p 



(4.8) 



We suppose the function ip(-) is finite on some interval 2 < p < b, where b = 
const < oo. It is evident that in the case b < oo ip(p) x ip(p). 

We introduce on the basis of the function ip(p) the new distance on the set T : 

1 4t,s) = \\g (t,r ] )-g°(s,r ] )\\G^ g ). (4.9) 

Further, put 

00) = = max sup log E (exp(e\g°{t,r]))) = 

£=±1 teT v ' 

max sup log / exp(e A g°(t,x)) u(dx). (4-10) 
£=±1 teT Jt 

If the function </>(•) is finite on some non-trivial interval A G (— Ao,Ao), where 
A = const > 0, we introduce on the basis of the function 

<f>g(\) = sup n(f>g(\/y/n) 

n=l,2,... 

the new distance on the set T : 

7*(M) = \\9%V) -9°(s,r,)\\B(<f, a ). (4.11) 



Obviously, 



7*(*, s) < Crfit, s), l4> (t, s) < C 2 P(t, s). (4.12) 



We denote for arbitrary separable numerical bounded with probability one ran- 
dom field ((t) = ((t,u); t G T, where (Q,B,P) is probabilistic space, 



P c (u)=p(sup|C(*)l>«), (4.13) 
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and correspondingly 



P (n) (w) := P ^Vnsup - I(t)\ > u^j 



We introduce also the centered separable Gaussian random process, more exactly, 
random field X(t) = X g (t), t G T with a following covariation function 

Z(t, s) = Z g (t, s) = cav(X(t), X(s)) = EX(t)X(s) = 

f g(t,x)g(s,x) u(dx)-I(t)I(s). (4.14) 

Theorem 4.1. Exponential level. 

Suppose the following condition is satisfied: 

C%-,(H(T,^e))de<oo. (4.15) 

J u 

Then the centered Gaussian field X(t) is continuous a.e. relative the distance 7^ 
and 

lim P (n) (w) = P x (u). (4.16) 



Theorem 4.2. Power level. 

Suppose the following integral converges: 

j\^{H{T^e))de< 00. (4.17) 

Then the centered Gaussian field X(t) is continuous a.e. relative the distance 7^ 
and there holds 

lim P (n) («) = P x (u). (4.19) 

n— >oo 

Remark 4.1. The exact asymptotic for the probability P^(u) as u — > 00 is 
obtained in [36], p. 19, 88, 106, 114, 180: 

P x (u) ~ C(X, T) M*" 1 exp (- M 2 /(2^ )) > 
where C(X,T) = const G (0, 00), 

a 2 + = max Z g (t, t) = max ^ g 2 (t, x)u(dx) — I 2 (t) , 

the value k = const dependent on the geometrical characteristic of the set 

T Q = {t,s:Z g (t,s)e[o 3 + /2,o 3 + ]}. 
A non-asymptotical estimation of Px(u) for u > 2a + of a view 
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Px(u) < C + (X,T) u^exp (-u 2 /(2a 2 + )) , 
is obtained, e.g., in [30], chapter 4, section 4.9. 

Remark 4.2. It is important for the practical using of offered method to calcu- 
late the main parameters o\ and k. It may be implemented by the following method, 
again by means of Monte-Carlo method. 

Let T = {t m } C T be some finite net on the whole set T. The consistent estima- 
tion cPl of a value o\ has a view: 



<7 + max Z g (t m , t m ) ~ max 



n _1 E^.6)-^ 

i=l 



Analogously may be computed the distances d^, d^(t,s) and D(t,s); for example, 
d(ti,t m ) w d(ti.t s ), where 



The consistent estimation of a value ||?y||5(0) based on the independent sample 
{r^j}, i — 1, 2, . . . , n is described, e.g., in the monograph [30], p. 291-294. 

This facts may be used by building of confidence region in the uniform norm 
for the integral I(t). Namely, let 5 be a "small" number, for example, 5 = 0.05 or 
5 = 0.01 etc. 

The value 1 — 5 may be interpreted as a reliability of confidence region. 
We define the value u(8) as a maximal solution of an equation 

Px(u(5))=5, 

or asymptotically equivalently, the maximal positive solution of an equation 

C(X,T) u(5) K - 1 exp (-u(5) 2 /(2al)) = 5. 

The asymptotical as 5 — > 0+ confidence interval in the uniform norm with reliability 
(approximately) 1 — 5 for /(•) has a view 



sup -/„(*) | 



(4.20) 



Proof of both theorems 4. Land 4.2. 

1. From the classical Central Limit Theorem (CLT) for the independent identi- 
cally distributed centered random vectors follows that the finite-dimensional distri- 
butions of the random fields 



X n (t) = VH(I n (t)-I(t)) 



(4.21) 



converge in distribution as n — > oo to the finite-dimensional distributions of Gaussian 
field X(t). It remains to prove the weak compactness of the set of (probabilistic) 
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measures in the Banach space of continuous functions C(T, 7^) or correspondingly 
C(T, 7^,) induced by the random fields X n (-). 

2. We use further the result of Lemma 3.2. Namely, we put A = 1,2, .. . and 
consider the differences 

[I n (t) - I(t)} - [I n (s) - I(s)] = n-^ ±[(g(t, fc) - 1(f)) " {9(8, fc) - /(a))]. 

i=i 

We conclude using the inequality (3.9) (and further inequality (3.8) 

\\[I n (t) -I(t)} - [I n (s)-I(s)]\\G(f g ) < C d^ g (t,s). (4.22) 

Since the exact value of constant C is not essential, we get to the assertion of theorem 
4.1. 

The proposition of theorem 4.2 provided analogously; instead the d^(-,-) distance 
we will use the metric •). 

5 Confidence region for solution of integral equa- 
tions 

Let us return to the source integral equation (1.1). We retain the notations of the 
sections 1 and 2: n,e,N = N(e),9(m),n(m) = 9(m) n,R(x),d(t,s),yf\t),y( N \t) 
etc. 

Another notations: 

Mp) = \R(-)\ p = 

and suppose ipn(-) G ^, i.e. ipn(b) < 00 for some b > 2; 
v^ m = inf 

2/,S/€(0,l),i/>fl(l/j/)<oo 

Vi/, = V^N = Vi(>,N(e), 

1(e) = / v^ N (H(T, d,x)) dx; (5.3) 

v 

Z m (t,s)= / J R'(t,a;i)ir(s,Xi) J R' 2 (xi,a;2)^ 2 (x2,X3) . . .i^ 2 (x m _i,x m )- 

/ 2 (x m ) n(dx 1 )fi(dx 2 ) ■ ■ ■ fJ,(dx m ), (5.4) 



\R(x)\ p n(dx) 



(5.1) 



xy + log 



V#(i/y) 

y-log(l + l/y) 



(5.2) 
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N 



Z(t,s) = ]T Z m (t,s)/9(m). (5.5) 

m=2 

Let X(t) = X e (t) be a separable centered Gaussian field with the covariation 
function Z(t,s) : EX(t) = 0, EX(t)X(s) = Z(t,s). 

A. We consider first of all the power level for integral equation. 

Theorem 5.1. a. Assume in addition that for some e G (0, 1) /(e) < oo. Then 
for such the value e the gaussian random field X(t) is D(-, •) continuous a.e. and 

hm P (v^max|yf \t) -y^ N \t)\ > «) = P (max|X(*)| > «) . (5.6) 

Theorem 5.1.b. Assume in addition that for arbitrary e G (0,1) /(e) < oo. 
Then for all the values e the gaussian random field X(t) is D(-, •) continuous a.e. 
and the proposition (5.6) holds. 

It is enough to prove only theorem 5.1. a. 

Proof consists on the using of theorem 4.1 to the each summand S™^[f](t). 
We need to prove only as before the weak compactness of the sequence of random 
fields 



Indeed, let \\f\\C(T) = 1 and let m = 1, 2, . . . , recall that Ve G (0, 1) N 
N(e) < oo. We have: 

\K(t,si)K(s 1 ,s 2 )K(s 2 ,s 3 ) . . .K(s m -i,s m )f(s m )\ < R(s 1 )R(s 2 ) ■ ..R(s m ), 

| [K(t, si) - K(s, si)] • K(s u s 2 )K(s 2 , s 3 ) . . . K(s m - U s m )f(s m )\ < 

d(t,s)-R(s 1 )R(s 2 )...R(s m ). 
By our condition, /?(•) G G(^p R ) or equally 

J T \R{ S )\^{ds)<r R {p)- 

/ \K(t,s 1 )K(s 1 ,s 2 )K(s 2 ,s 3 ) . . .K(s m _i,s m )f(s m )\ p ]J n(ds k ) < 
J1 k=i 

m m m ,. 

n R p ( s k) n v{ds k ) = n / ir^M^) < r tp w 
k=i k=i k=i jT 



Hence 



and analogously 
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„ m 

/ | [K(t, Sl ) - K(s, Sl )} ■ K( Sl , s 2 )K(s 2 , s 3 ) . . . K{s m ^ s m )f(s m )\P J[ »(ds k ) < 
jTm k=i 

dP{t,s)^{p). 
Therefore, the random process ((t) = ( m (t) = 

K(t, zjKfo, &m 2 , . . . K(t m -i,uf(u 

belongs to the space G(ip m ) uniformly on t G T : 

sup||C m (t)||G(^ m )<l 
teT 

and 

\\ut)-Us)\\G(r i )<d(t,s). 

The application of theorem 4.1 completes the proof of theorem 5.1. 
B. Integral equations. Exponential level. 
Let us define the function 0(at)(-) from the set $ : 



4>{N)(p) = 



P 



N 



-1 



,p>2, (5.7) 



7r(A) = 7Tat(A) = sup [m<f>w)(^/ Vm)]- (5.8) 

m=l,2,... 

and suppose the finiteness of such a functions for some values p > 2. 
Denote ( 

J=J{e) = [ xAH(T,d^x))dx. 
Jo 

Theorem 5.2. Assume in addition that for any e e (0, 1) J — J(e) < oo. Then for 

such the value e the gaussian random field X{t) is D(-, •) continuous a.e. and the 
proposition (5.6) holds. 

Proof is at the same as in theorem 5.1a and may be omitted. 



6 Non-asymptotical approach 

We evaluate in this section the non-asymptotical probabilities for deviations 

P(w) = supP (n) (w) = supP(v^sup|/ n (t) -I(t)\ > u) (6.1) 

n>l n>\ teT 



and correspondingly 



Q(u) = supQ (n) (w), 

n>l 
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where 

Q^(u) = P (vWx \ y f \t) - y W(t)\ > «) = P (max \Y n (t)\ > u) , (6.2) 

Fn^^^rW-y^W)- (6-3) 

A. Multiple integrals. Exponential level. 
Assume for some function = 0(A) G $ 

sup | \g(t, V ) - I(t)\\B(<f>) = f < oo. (6.4) 

The condition (6.4) is satisfied, e.g., for the natural choice of the function 0(A) = 
0o(A) and = 1. 
Recall that 

d*(t, s) = 1 1 [g(t, V ) - I(t)} - [g(s, rj) - I(s)} \ \B (</>). 
Theorem 6.1a. Suppose 

J ^) : =j X^(H(T,d (l> ,x)dx <oo. (6.5) 



Then for the values u > 2y J(0) 

P(a^u) < 2 exp (-0* (u - ^2 J(0) u)J . (6.6) 
Proof used the lemma 3.1. We get using the definition of the function 0(-) : 
sup ||v^(/ n (t)-/(t))||£?(0)<^; 

n=l,2,... 

SUp 1 1 [(/„(*) "/(*)) - (InW-Hs))] ||5(0) < d+frs). 
n=l,2,... 

This completes the proof of theorem 6.1. 
B. Multiple integrals. Power level. 

Theorem 6.1b. Suppose for some function ip — ip(p) e ^(2, b), b > 2 

sup | \g(t,ri) - I(t)\\G(i/j) = f a* < oo. (6.7) 

teT 

The condition (6.7) is satisfied, e.g., for the natural choice of the function ip(p) = 
ip (p) and cty = I. 
Recall that 

d f (t,s) = \Mt,r,)-I(t)] - [g(s,r ] )-I(s)]\\GW. 

We define 
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Z($) := + 9 J (\og(2N(T, d f , x))) dx < oo. (6.8) 

If Z(ip) < oo, then for the values u > 2Z(ip) 

P(«) < exp {-w^u/Z^))) . (6.9) 

Proof is at the same as the proof of theorem 6.1; it used the lemma 3.2 instead 
lemma 4.1 and the definition of the function ip(-). 
For instance, if 

</>o (A) ~ A r , A > 1, r = const > 1, 
and J ((f)) < oo, then for u > 1 

P(u)<exp (-d(r, J(V>)K), 

f== m ;; (2 ; r) , (6.io) 

mm(2,r)-l v ' 

C. Integral equations. Power level. 

Analogously to the theorem 6.1a may be proved the following two results. 
Theorem 6. 2. a. Assume that for some e G (0, 1) /(e) < oo. Denote 

We deduce: 

Q(u) < exp (-«,; w (CfV(l + He)))) ■ (6-11) 
If for example X = R d , 

fi{x : R(x) > u} < exp (-C 2 u l,p ) ,u > 1, 
and T is bounded open subset R d , 

d^(t, s) < C 3 (min (J log \t - s\ |~ 7 , l)) , 

7 = const >^ = / /3A^(e) + l, 
then Z(ip) < oo and following 

Q(u) < exp (-C 4 (C 2 ,C 3 ,7,£,d) u 1 ^) ,« > 1- (6-12) 



Remark 6.1. We accept that the case /3 = is equivalent the boundedness of 
the function R(-) : 

/3 = <^ vraisup i?(x) < oo. 
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D. Integral equations. Exponential level. 
We define the function 0(at)(-) from the set $ : 

-l 

,p>2, (6.13) 
7r(A) = 7rjv(A) = sup [m(j)f N \(X/y/m)]. (6.14) 

m=l,2,... 

and suppose the finiteness of such a functions for some values p > 2. 
Denote 

Jo 

and suppose J < oo. 
Proposition: 

Q(«) < exp (— 7r*(u - v 7 ^)) ,u>2J. 
As an application: solving the equation 

F(«p(5)) = 

or correspondingly 

QK>(<5)) = 5 

relative the variable u = u P (S) or u = uq{8) where as before 1 — 5, 5 = 0.05, 0.01 
etc., is the reliability of non- asymptotic confidence region in the uniform norm, we 
conclude that with probability at least 1 — 5 

I(t)\<u P (5)/V^, (6.12) 
y^(t)\<u Q (5)/V^. (6.13) 



<t>{N){p) 



P 



N 



sup \I n (t) - 
teT 



7 Examples 

We suppose in this section that T is bounded closed domain in the space R d with 
positive Lebesgue measure n(D) = J D dx. Denote as ordinary by \t — s\ the Euclidean 
distance between a two points t, s;t,s e T. 
We assume again r(S) < 1, r(U) < 1. 

Example 7.1. Multiple parametric integral. "Power" level. Recall that 

g(t,x) v(dx); I n {t)=n- 1 Y,9(t,r]®), 
1 i=i 

and Law (77(2)) = ia 
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Assume that for some a G (0.1] 

\g(t,x)-g(s,x)\<\t-s\ a Q(x) (7.1) 

where for some 5 > 

/ Q d ' a+5 {x) u(dx) < oo & Q(-) e L d/a+s . (7.2) 

Since 

N(T, \t-s\ a ,e) -e- d/a , e -)• 0+ 
we conclude that for the value p , where 

— = d/a + 5, 
Po 

the Pizier's condition of theorem 4.1 is satisfied. 

Example 7.2. Multiple parametric integral. "Exponential" level. Here 
we suppose 



\g(t,x) -g(s,x)\ < max[| log |t - s\\ 7 , 1] Q(x), (7.3) 



where 



v{x : Q(x) >u}<d exp [-C 2 u 1/I3 ) , (7.4) 

and 

/3,7 = const, 7 > (3. (7.5) 

We conclude that the condition of theorem 4.1 is satisfied. 

Note that the condition (7.3) is equivalent to the following inequality: 



P>1 

or equally 



sup \Q\ p /p p < oo, 
p>i 



Q(-)GG(^), ^(p) =V- 
The condition of theorem 4.1 is satisfied also when 

sup \Q\p/ p 13 (hgp) 132 < oo, f3 2 = const > 0. 
p>i L J 

Example 7.3. Integral equation. "Power" level. 

A. Theorem 5.1a. 

Let again T is bounded open subset of the space R d and assume as before 

fi{x : R(x) > u] < exp (-du 1 ^) , u > 1, (7.6) 
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m*, s ) < c * ( min (i lo s i* - s i r 7 . !)) > ( 7 - 7 ) 

7 = const > /3N(e) + 1. (7.8) 
Then all the conditions of theorem 5.1a are satisfied. In particular, Z(ip,e) < oo. 
B. Theorem 5.1b. 

Let the condition (7.6) be satisfied. Suppose also (instead conditions (7.7) and 
(7.8) ) 

d(t, s) < C 4 \t - s\ a , a = const £ (0, 1]. 

Then all the conditions of theorem 5.1b are satisfied. 
In particular, Ve £ (0, 1) Z(tp, e) < oo. 



8 Derivative computation 

Let us return to the source equation (1.1). We consider in this section the case 
T = [0, 1] with the classical Lebesgue measure /i(dx) = dx. Suppose / be continuous 
differentiable: // £ C[0, 1] and that there exists a continuous relative the variable t 
function 

We obtain denoting Y(t) = yl (t) after the differentiation of equation (1.1) 

fi 
io 

We do not entrust the contraction condition on the kernel V = V(t,s). 



Y(t) = f'{t) + / " V(t, s)y(s)ds d = f f~(t) + V[y](t). (8.1 



The solution Y(t) may be written as follows. 

oo 

Y(t) = f(t) + J2 °m(t), 
m=l 

where 

vi{t) = / V(t,x)f(x)dx, m = 2,3, . . . =>• a m (t) = 
Jo 

f ds f V(t,s)K(s,s 1 )K(s 1 ,s 2 ) ■ ■ .K(s m -i,s m )f(s m ) ds ds 1 ds 2 ■ ■ -ds m . (8.2) 

JO JT m 

Let the number e, e £ (0, 1) be a given. We retain here the notations of the 
second section: N = N(e), n(m) = 9(m) ■ N(e) etc. 

We choose as the deterministic approximation Y^ N \t) for Y(t) as before the 
expression 
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N(e) 

Y( N \t) = f~(t)+J2a m (t). (8.3) 

m=l 

The accuracy of the approximation (8.3) ("bias") in the uniform norm may be 
estimated as follows. 

sup \YW(t)-Y(t)\<\\V\\- sup \yW(t)-y(t)\<\\V\\-e. (8.4) 
te[o,i] te[o,i] 

Recall that 

||V"|| = sup / \V(t, x)\dx < oo. 

*6[0,1] JO 

Further, we offer for the Y^ N \t) calculation the following Monte-Carlo approxima- 
tion. Let us consider separately the expression for a m (t), m — 2, 3, . . . , N. Let the 
random variable ( be uniform distributed in the set [0, 1] and let the random vector 

— * 

£ = {£i?£2, • • • ,£m} be uniform distributed in the multidimensional cube [0, l] m ; let 
also {C (j) }, {fW} = {fi, &...,&}, 3 = l,2,...,n(m) be independent copies of 
the vector {C, £}• 
Note that 

a m (t) = EV(t, ()K((, 6) • • • K(£ m -i, U)f(U)- (8.5) 

We can offer therefore the following Monte-Carlo approximation a m (t) for a m (t) : 

n(m) 



1 



X 



/Y(c^,d J) )^(d J) ,#) • . .tftfSL.eff) (8.6) 

and correspondingly the following approximation Y^ N \t) for Y^ N \t) : 

N(e) 

Yt N) (t) = f~(t)+j:a m (t). (8.7) 

m=l 

Here n denotes as in the section 2 the common quantity of elapsed random variables. 
Note that 

Var m = Var[a m (t)] < — — ds ds 1 ds 2 . . . ds m V 2 (t, s)f 2 (s m )x 

n(m) Jo JT m 

K\s,s l )K\s 1 ,s 2 )...K\s rn ^s m )<^—\\V^\\ ■ \\U m \\. (8.8) 

n(m) 

It is easy to calculate as in the second section that 
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N(e) 

supVar[F n (7V) (t)] = J2 supVar m < C/n. (8.9) 
tax m=l ter 

The inequality (8.9) show us that the speed of convergence Y r [ N \t) to Y^ N \t) is 
equal to 1/ ^Jn in each fixed point to £ T. 

In order to establish this result in the uniform norm for the sequence of random 
fields Y^ N \t), we need to introduce some new notations and conditions. 

\V(t, x) — V(s x)\ 

Ry(x) = vraisup \V(t, x)\, d v (t, s) = vraisup — — , 

ter xeT Ry{x) 

ip v (p) = \R v (-)\p, 

so that 

\V(t,x)\<R v (x), \V(t,x) -V(s,x)\ < d v (t,s) R v (x), (8.10) 

3by > 2, Vp < by 4>v(p) < oo; 
buy = min(6, by) = const > 2; 

V^t, s) = j T V(t,z)V(s,z)dz; 
V m (t, s)= dz V(t, z)V{s, z)K 2 {z, x^K^xuXi) . . . K 2 (x 

m—\i -Km) 

JO JT m 



X 



f 2 {x m ) dx x dx 2 ■ ■ ■ dx m , m = 2, 3, . . . , N(e); (8.11) 



AT(e) 1 



Zy(t,s) = Zy iN (t,s) = Fl^f 4 '^ 



m= 
N 



®v(p) = Mp) ■ [Mp)Y , pe(2,M; 

Zy.= 1 + 9 C v* ev (H(T,d v ,x)) dx; (8.12) 
Jo 

Let Xy(t) = Xy tt (t) be a separable centered Gaussian field with the covariation 
function Z v (t,s) : EX v (t) = 0, EX v (t)X v (s) = Z v (t,s). 

Theorem 8.1. (Power level.) Let Zy < oo. Then the limiting Gaussian random 
field Xy(-) is d v continuous a.e. and 



lim P(v^sup \Y r [ N \t) - Y {N \t)\ >u)= P(sup \X v (t)\ > u), u> 0; (8.13) 

n-HX, t( - T t( - T 
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supP( v / ^sup \Y^ N) (t) - Y^ N \t)\ >u)< exp (-w* e (u/Z v )) , u > 2Z V . (8.14) 



Proof is at the same as the proof of theorems 5.1a and 5.1b. 

Namely, let us introduce the following random processes (fields) 

Cm,v(t) = V(t, ()K((, tfKfa, 6) • • • K(U-u U)f(U~ 



V-S m [f](t), m = l,2 



N 



so that E( m y(t) = 0, and its independent copies 



(8.15) 



V-S m [f}(t), m = l,2,...,JV; 



-. n(m) 

1 ^ -0') 



E Cm 



(8.16) 



(8.17) 



>n(m) j=i 

#w=c#,vw, s(«)(t)=H?(t). 

It is sufficient to consider only the case m = N. We need to prove the tightness the 
random fields S^(t). 

We have assuming without loss of generality sup t \f(t)\ = 1 : 

\(m,v(t)\ < Rv(0 R(£i)Rfa)...R(U), 

\(m,v(-)\ P <Mp) 

and using the Rosenthal's inequality 



\Z {n) (t)\ p <C^p MP) «p)/logP- 

Analogously 

|S<">(f) - E^(s)\ p < d v (t, s) ■ C^p MP) «p)/logP, 
We conclude after summing over m : 



sup sup 

n t 



V^(Y^ N \t) -Y^(t))]\\G(e v ) <1, 



(8.18) 



(8.19) 



(8.20) 



sup 

n 



^(t) " y (N \t))] - [V^M N \s) - YW(s))] \\G(Qy) < 

d v (t,s). (8.21) 
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The assertion of theorem 8.1. follows now from lemma 3.2 and theorem 3.2. 
Example 8.1. If 

fi{x : R(x) > u} < exp , (3 = const > 0, 

fi{x : Rv{x) > u} < exp (— Cu 1 ^^ , u = const > 0, 
d v (t,s) < C 

and 

7 > N(e)p + uj + l, 
then the conditions of theorem 8.1 are satisfied. As a consequence: 

sup P(^ sup \Y^ N \t) - Y^ N \t)\ >u)< exp (-C((3, ^u 1 ^ 1 ') , u > 2. 

(8.22) 

In the case when 

d v (t, s) < C\t — s\ a , a = const > 0, 

then the conditions of theorem 8.1 are satisfied for arbitrary values e£ (0,1). 
Remark 8.1. CLT in the space C^T]. 

We used in this section in fact the Central Limit Theorem in the space of con- 
tinuous differentiable function C^T]. 

Remark 8.2. Some generalizations. Let A = A t be any linear operator, not 

necessary to be bounded, defined on some (dense or not) subspace of the space C(T), 
for example, differential operator A = d/dt, partial differential operator, Laplace's 
operator etc. 

We write formally 

A t y = A t f{t) + / A t K(t, s)y(s)fjL(ds). 
Jt 

Using at the same considerations, we might obtain the CLT in the space Ca[T] 
consisting on the continuous functions g = g(t) with continuous Ag(t) equipped by 
the "energy" norm 

= sup \g(t)\ + sup \Ag(t)\. 

teT teT 

We can conclude as a consequence that the rate of convergence of the random 
approximation for the Monte-Carlo approximation 



mm 



M log |t — s| | 7 j , 1 ,7 = const, 



41 



JV(6) 

Yi N \t) :=Af(t) + J2Aa m (t) 

m=l 

in the space Ca [T] is equal to 1/ y/n and the bias is less than C ■ e : 

HAy-yf^lU^e-sup / \A t K(t,s)\ds. (8.23) 

We consider now the exponential level for the derivative computation. Recall (see 
section 3) that the estimations through B(<f)) spaces (exponential level) have advan- 
tage in comparison to the estimations using the G(ip) technique if the correspondent 
</>(•)— function there exists. 

Define a function 

(8.24) 

if there exists for some interval p G [1,6], b = const > 2, and suppose (py(-) G 
and introduce the correspondent distance 

<W(M) = \\(N,v(t) -CN,v(s)\\B(<j>v). (8.25) 
Moreover, assume that the following integral converges: 

K v := f ' x^{H{T,d^,x)) dx < oo. (8.26) 

J 

Theorem 8.2. (Exponential level). 

Let Ky < oo. Then the limiting Gaussian random field Xy(-) is d$ v continuous 
a.e. and 



(hip) = 



V 



<3>v(p) 



lim P(v^sup \Y^ N) {t) - Y iN \t)\ > u) = P(sup \X v (t)\ > u), u> 0; (8.27) 
supP( v / ^sup \Y r [ N \t) - Y {N) (t)\ >u)< exp (-<h(u - \/2K v u)) , u > 2K V . 

n teT V / 

(8.28) 

Proof is at the same as the proof of theorems 8.1 and may be omitted 

Remark 8.3 As we know in the remark 4.1, the tail functions for the maximum 
distributions in (5.6) and (8.27), i.e. probabilities 

P(max|X(t)| > u), P(max \X v (t)\ > u) 
have the exact asymptotic of a view, e.g., 

P(max|X(t)| > u) ~ C(X,T) u K ~ l exp (-u 2 /{2a 2 + )) , 
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where a\ = max t Var[X(t)] etc. 

The consistent estimations of the parameters k, <j\ and C(X,T) may be the 
implemented as in the remark 4.2, where we can computed multiple integrals by 
Monte-Carlo approximation. 

Obviously, these parameters may be estimated through the analytical expression 
for the functions f,g,K. 

Example 8.2. The conditions of theorem 8.2 are satisfied if for example both 
the r.v. R(x), R v (x) are essentially bounded, X is bounded subset of the space R d 
and 

<W(*> s ) ^ °5 ( min (l lo s I* _ s l I -7 ' x )) ' 7 = const > 2 ( 8 - 29 ) 

or moreover in the case when (instead the condition (8.29)) 

<W(*> s ) <Ce\t- s\ a , 3a = const e (0, 1]. (8.30) 



9 Concluding remarks 

A. Another method. 

Let us consider the following equivalent modification of source equation (1.1): 

y x (t) = f(t) + \J t K(t, s) y(s) fi(ds) = f(t) + XS[y](t), (9.1) 
where A = const G (0, 1) and as before 

Pl = r {S) < 1, p = r(U) = r (S (2) ) < 1 

and suppose that the function /(•) and the kernel K(-, ■) satisfies all the conditions 
of the section 1. 

The solution y\ = y\(-) may be written as follows: 

oo 

yx = f+T, xmsm [f}- ( 9 - 2 ) 

m=l 

Let us introduce the so-called geometrical distributed integer random variables 

t : 

P(r = m) = (1 - A)A m , m = 0, 1, 2, . . . ; 

then 

(l-X)y x = ES T [f}. (9.3) 

Let M — 2, 3, ... be arbitrary integer number and let Tj, j = 1, 2, . . . , N be inde- 
pendent copies of r. It may be offered as an consistent approximation for (1 — \)y x 
the following expression: 
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M 

{l-\)y x , n {t) = M' l Y J S^\ 
where the value has a following Monte-Carlo approximation: 

n(j) 

S rU) „ gr U ) ^ (n( j))-l £ ^O-))^^). 

i=l 



So, 



M n(j) 

;i - \)yx, n = m- 1 2 s^Wj))- 1 E ^ (T(j)) [/](SVw)- ( 9 - 4 ) 
i=l 1=1 



We denote also 



n:= (n(l),n(2),...,n(M)); 



n is any M— tuple of integer positive numbers. 
We obtain after some calculations: 



1 ™ 1 



<KM,n):=Var(y Ain )x— £ 



Note as before that the common amount of used T— valued random variables is 
equal to 

M 

A(M, n) :=M.-£jn(j). 
We conclude solving the following constrained extremal problem: 

(f>(M, h) — > min /A(M, h) = n, 

that the minimal value of the function <f)(M,n) under condition A(M,n) = n is 
asymptotically as n — > oo equivalent to 

min 0(M, n) /[A(M, n) — n] x n - 1/2 . (9.5) 

Thus, the optimal speed of convergence y\ )Tl to the solution y is asymptotical equal 
to n" 1//4 , in contradiction to the first offered method. 

Another sufficient conditions for CLT in the space of continuous functions C (T) 
see, e.g. in the works [15], [20], [38], [49]. 

We adapt here only the results belonging to Marcus M.B. and Jain N.C. for the 
integral I(t) = f x g(t,x)v(dx) calculation. Namely, if 

\g(t,x)-g(s,x)\<M(x)-p(t,s), (9.6) 
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(condition of factorization), 

f M 2 (x)u(dx) < oo, (9.7) 

J X 

(moment condition), where p(t, s) is some metric on the set T for which 

f 1 H 1/2 (T,p,z) dz < oo, (9.8) 
Jo 

(entropy condition), then the sequence of random processes (fields) y/n(I n (t) — /(£)) 
satisfies the Central Limit Theorem in the space C(T,p). 
Analogously, if 

( exp A ([g(t, x) - g(s, x)} - [I(t) - I(s)}) v{dx) < 

J X 

exp(0.5A 2 A 2 r 2 (t, s)), A = const > 0, (9.9) 

or equally 

|| ([g(t,x)-g(s,x)} - [I(t) - I(s)}) || sub < Ar(t, s) 
(subgaussian condition), for some (semi-) distance r(-, •) for which 

f 1 H 1/2 (T,r,x)dx < oo, (9.10) 
Jo 

(entropy condition), then also the sequence of random processes (fields) y/n(I n (t) — 
I(t)) satisfies the Central Limit Theorem in the space C(T, t). 



10 Necessity of CLT conditions 

We will confer in this section the necessity of some conditions for the Central Limit 
Theorem in the space of continuous functions C(T,p). 

A. Condition of factorization. 

Let £(t) be continuous with probability one random field relative some distance. 
Then there exist a non-random continuous distance p = p(t, s) and a random variable 
M = M(u) for which 

\((t)-((s)\<M(u)-p(t : s) : (10.1) 

see [47], [5]. 

B. Moment condition. 

Assume in addition the field ((t) satisfies the Central Limit Theorem in the space 
of continuous functions. Then £(•) has the weak second moment. Following, see [47], 

EM 2 < oo. (10.2) 

C. Entropy condition. 
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Let the conditions (9.6) and (9.7) for the field ( = ((t) be satisfied; assume 
also the field ((t) satisfies the Central Limit Theorem in the space of continuous 
functions. 

Suppose in addition T = [0, 2ir] d , d = 1, 2, . . . and that the distance p(t, s) from 
the inequality (10.1) is translation invariant: 

p(t,s) = p(\t-s\), (10.3) 

where t ± s — t ± s(mod27r). 

Other assumption. Let the random field ((t) satisfies the Central Limit Theorem 
in the space of continuous functions. 

Denote by CooW the separable centered Gaussian field with at the same covari- 
ation function as r](t) : 

ECocWCoc(s) = EC(f)C(*), (10.4) 

then 

/UM) :=£!/%(*) -Cco(^ * 

E 1/2 [CW-C(s)] 2 =:rf C (t, S ). (10.5) 
It follows from the condition (10.2) that 

d c (t,s)<Cp(\t-s\); 
assume in addition the conversely: 

p(|t-s|)<CH(M). (10.6) 
As long as the distance p<x,(t, s) is linear equivalent to the distance p(\t — s\) : 

d c (t,s) x p(\t - s\),t,s G T, 

the convergence of entropy integral (9.8) follows immediately from the famous result 
of X.Fernique [11]. 

D. Subgaussian condition. 

Let the condition (9.9) be satisfied. Suppose alike in the last pilcrow in addition 
T = [0, 27r] d , d — 1,2,... and that the distance r(t, s) from the inequality (9.9) is 
translation invariant: 

r(t,s)=r(\t-s\), 

where t ± s — t ± s(mod27r). 

Then the (weak) limiting Gaussian field v = v(t) in the CLT for the space of 
continuous functions satisfies at the same condition (9.9). 

It is evident that 

\u(t)-u(s)\ 2 <T(\t-s\); 
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assume conversely, i.e. that 



r(\t - s\) < C\u(t) - u(s)\ 2 ; 

then the convergence of entropy integral (9.10) follows again from the result of 
X.Fernique [11]. 
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